Black holes and stars in Horava-Lifshitz theory with projectability condition 
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We systematically study spherically symmetric static spacetimes filled with a fluid in the Horava- 
Lifshitz theory of gravity with the projectability condition, but without the detailed balance. We 
establish that when the spacetime is spatially Ricci flat the unique vacuum solution is the de Sitter 
Schwarzshcild solution, while when the spacetime has a nonzero constant curvature, there exist two 
different vacuum solutions; one is an (Einstein) static universe, and the other is a new spacetime. 
This latter spacetime is maximally symmetric and not flat. We find all the perfect fluid solutions 
for such spacetimes, in addition to a class of anisotropic fluid solutions of the spatially Ricci flat 
spacetimes. To construct spacetimes that represent stars, we investigate junction conditions across 
the surfaces of stars and obtain the general matching conditions with or without the presence of 
infinitely thin shells. It is remarkable that, in contrast to general relativity, the radial pressure of a 
star does not necessarily vanish on its surface even without the presence of a thin shell, due to the 
presence of high order derivative terms. Applying the junction conditions to our explicit solutions, 
we show that it is possible to match smoothly these solutions (all with nonzero radial pressures) to 
vacuum spacetimes without the presence of thin matter shells on the surfaces of stars. 

PACS numbers: 04.60.-m; 98.80.Cq; 98.80.-k; 98.80.Bp 
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I. INTRODUCTION 

There has been considerable interest recently on a non- 
relativistic theory of gravity proposed by Horava [l| . In- 
spired by the theory of a Lifshitz scalar relevant in 
condensed matter systems, Horava postulated a nonrela- 
tivistic anisotropic scaling symmetry of space and time. 
This allows for the addition of higher spatial derivative 
terms in the action without their time derivative counter- 
parts, rendering the theory power counting renormaliz- 
able. The theory is supposed to flow dynamically from a 
scale invariant theory in the ultraviolet (UV) to General 
Relativity (GR) in the infrared (IR), thus restoring the 
diffeomorphism symmetry at low energies. 

In the class of Horava-Lifshitz (HL) theories studied 
during the past year, there have been many variations 
regarding both the potential that is used and whether 
the lapse function is allowed to depend on the spatial co- 
ordinates or not. In the initial proposal by Horava [lj], 
a simplifying assumption was made to reduce the num- 
ber of terms in the potential. The potential was required 
to be derived by a "superpotential," thus giving is it a 
form dubbed as detailed balance. However, it was soon 
realized that the breaking of detailed balance is neces- 
sary in order to obtain Minkowski vacua at low energies 
0, 0] . This breaking can be done minimally by adding a 
soft breaking term in the action, or can be done in more 
generality by adding all possible terms. In the present 
paper, we will adopt the potential of [3j, where all oper- 
ators that conserve parity are included in the potential. 
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Since the detailed balance potential has operators that 
break parity, there will be no limit of the dynamics of 
the potential constructed in Q, unless the contributions 
of these operators vanish. It can be shown that when 
the spacetime is static and spherically symmetric, these 
operators indeed have zero contributions Q. A second 
assumption that distinguishes HL theories, is whether 
one assumes that the laspe function A is a function only 
of time (projectable case) or a spacetime function (non- 
projectable case). In the original proposal, projectabil- 
ity was assumed, but since the Hamiltonian constraint is 
not local (it becomes an integral constraint Q), it was 
preferred in the literature to use the nonprojectable as- 
sumption. However, the latter seems to be inconsistent 
[6| as the Poisson brackets of the theory do not form 
a closed structure. For this reason, in the present pa- 
per we choose to work with the projectable assumption. 
As is expected, the breaking of diffeomorphism invari- 
ance results in the existence of an extra (scalar) mode 
in the spectrum, which is absent in GR. This mode in- 
troduces pathologies for almost all variations of the the- 
ory, giving rise to sometimes instabilities (classical or 
quantum), and always to strong coupling. In particu- 
lar, in the projectable case this mode exists already in 
Minkowski backgrounds and is either classically or quan- 
tum mechanically unstable @, 01- I n the nonprojectable 
case, this mode exists for time-dependent and spatially 
inhomogeneous backgrounds and is classically unstable 
8] . If one includes spatial gradients of the lapse function 
in the nonprojectable case [9(, the mode can be rendered 
stable. The cases described above that are classically un- 
stable can be rendered stable by higher order interactions 
with the price of introducing strong coupling in the the- 
ory Jj| [T3| • Strong couplingalso exists in the limit that 
the theory approaches GR [9] , which brings the effective 
perturbative scale lower than the GR cutoff [ill ]. 
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To understand the theory further, other aspects have 
also been studied [12]. In particular, solutions of this the- 
ory have been extensively investigated in the directions 
of cosmology as well as black holes. Isotropic cosmolog- 
ical solutions were studied in [l3j and revealed that the 
new terms added in the theory typically modify the dy- 
namics for nonzero spatial curvature. It is interesting 
to note that in [5j it was advocated that the projectable 
theory has a dark matter component built-in, due to the 
nonlocality of the Hamiltonian constraint. 

The other set of solutions studied so far are the spher- 
ically symmetric ones [14j . In particular, Lu, Mei and 
Pope found all the vacuum solutions of the HL theory 
with the detailed balance condition 15]. Cai et al gen- 
eralized them to topological black holes with or without 
charges, and paid particular attention to their thermo- 
dynamics [l6| , while Colgain and Yavartanoo obtained a 
class of dyonic solutions with detailed balance and assum- 
ing that high order derivative terms in the potential of 
the massless vector field were absent [H (See also (3.). 
Adding a linear term of the three-spatial curvature into 
the action, so that the detailed balance condition was 
softly broken, Kehagias and Sfetsos constructed a class 
of vacuum solutions that is asymptotically flat [4] . Park 
generalized it to the case with a nonvanishing cosmolog- 
ical constant [19j . while Lee, Kim and Myung studied 
A1IS2 x S 2 solutions in such a deformed generalization 
of the HL theory [20|. Capasso and Polychronakos con- 
sidered the case with nontrivial lapse function and shift 
vector and found all the solutions with soft breaking of 
the detailed balance [2l[. Kiritsis and Kofinas further 
generalized the HL theory to include quadratic terms of 
the three-dimensional spatial curvature, and found the 
most general solutions 22], while Kiritsis himself [23[ 
considered static spherically symmetric solutions in a ver- 
sion proposed in [9j. Other kinds of solutions can be 
found in [12, [H| • It should be noted that all the solu- 
tions mentioned above were constructed without assum- 
ing the projectability condition. As we shall show in 
this paper, all such solutions can be written in a canon- 
ical ADM form that exhibits explicitly the projectabil- 
ity condition. In addition, static spherically symmetric 
spacetimes with the projectability condition were studied 
recently by Tang and Chen, and some solutions with or 
without charge were found [2J] . 

In this paper, we systematically study static spheri- 
cally symmetric solutions in the projectable case, with 
the full potential of @. We first present a self-contained 
introduction of the theory in Sec. II, and then in Sec. 
Ill we write down the equations of motion for projectable 
gauge that we will work with. In Sees. IV and V, we 
present all the vacuum solutions as well as solutions in 
the presence of a perfect fluid for two cases: the case 
that the spatial curvature vanishes (Sec. IV) and the 
case where the spatial curvature is a nonzero constant 
(Sec. V). In the former case we also present a class of 
solutions that represent an anisotropic fluid whose radial 
pressure is proportional to its tangential one. In Sec. VI, 



we first develop the general junction conditions across 
the surface of a star with or without the presence of an 
infinitely thin shell. Then, we join our vacuum solutions 
with the ones of a fluid by requiring that no thin shell 
be present on the surface of the star. Finally, in Sec. VI 
we present our main results. Three appendices, A, B and 
C, are also included. In Appendix A, the general expres- 
sions of the (F s )ij tensors appearing in the equations 
of motion for static spherically symmetric spacetimes are 
given, while in Appendix B static solutions in GR are 
studied in the canonical Arnowitt-Deser-Misner (ADM) 
form. In the latter, we show explicitly how one can 
bring any given static metric into an ADM form with 
the projectability condition. However, the kind of coor- 
dinate transformations needed to achieve this, are not al- 
lowed by the foliation-preserving diffeomorphisms of the 
HL theory, and the actions are generally not invariant. 
In Appendix C, we calculate the trace of the extrinsic 
curvature K — K\ for all the solutions found in Sees. 
IV and V, and study its singular behavior. 

Before proceeding, we would like to note that spheri- 
cally symmetric spacetimes in the framework of the HL 
theory with the projectability condition are also studied 
recently by Izumi and Mukohyama [25j, and found that, 
among other things, globally static and regular perfect 
fluid solutions do not exist. Our results presented in this 
paper do not contradict to it, since in this paper we have 
different assumptions. In particular, the perfect fluid to 
be considered here usually conducts heat al ong the radial 
direction [26[, while the one considered in [25j does not. 
This is another peculiar feature of the HL theory: Static 
stars in GR do not conduct heat. In addition, black holes 
in the HL theory mi ght not black, because of the different 
dispersion relations |22l [25| . 

II. HORAVA-LIFSHITZ GRAVITY WITHOUT 
DETAILED BALANCE 

We give a very brief introduction to HL gravity with- 
out detailed balance, but with the projectability condi- 
tion. (For further details, see [1,0, H3 ) The dynamical 
variables are N, N l and g^ (i, j = 1,2,3), in terms of 
which the metric takes the ADM form, 

ds 2 = ~N 2 dt 2 + g i:j (dx l + N l dt) (dx j + N j dt) . (2.1) 

The theory is invariant under the scalings 

t -> £ 3 t, x i -> lx i , 

N -> r 2 N, N l -> r 2 N\ gij -> g t3 . (2.2) 

It should be noted that there was a constant term c 2 in 
front of the lapse function N in the ADM metric used 
in |3j, so that N was rescalingas iV — > N, where c has 
dimensions [c] = [dx/dt]. In [3, H3| and this paper, we 
absorb it into N so that now N is scaling as that given 
above. The projectability condition requires a homoge- 
neous lapse function: 

N = N(t), N l =N l (t,x k ), n,, //,,(/../). (2.3) 
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The form of the metric is invariant under the foliation- 
preserving diffcomorphisms of the HL theory, 

t = t + X °(t), x i = x i + X i (t,x k ). (2.4) 

The total action consists of kinetic, potential and mat- 
ter parts, 

S = ( 2 J dtd?xN^/g (C K -C v + C 2 Cm) , (2.5) 

where g — det <7y , and 

C K = KijK* -(l-£)K 2 , 
C v = 2A - R + 4 (. g2 i? 2 + gsRijRV) 



c 2 



Q 4 



[g 7 RV 2 R + g 8 (V^fc) (V l i? jfc )] ■ (2.6) 



Here C 2 = l/167rG, the covariant derivatives and Ricci 
and Riemann terms are all constructed from the three- 
metric gij, while Kij is the extrinsic curvature, 

Kij = ^ {-gij + ViN, + VjNi) , (2.7) 

where Ni = gijN^ . The constants £, gj (I = 2, . . . 8) are 
coupling constants, and A is the cosmological constant. 
It should be noted that Horava included a cross term 
CijR^ , where Gy is the Cotton tensor. This term scales 
as £ 5 and explicitly violates parity. To restore parity, this 
term was excluded in Q. 

In the IR limit, all the high order curvature terms (with 
coefficients gi) drop out, and the total action reduces 
when £ = to the Einstein-Hilbert action. 

Variation with respect to the lapse function, N(t), 
yields the Hamiltonian constraint, 

j d 3 x^(C K + C v ) = 8^G J d 3 x4~g3\ (2.8) 



where 



J = 2 \ N 



.SC 



M 



SN 



C 



(2.9) 



Because of the projectability condition N = N(t), the 
Hamiltonian constraint takes a nonlocal integral form. If 
one relaxes projectability and allows N — N (t, x l ) , then 
the corresponding variation with respect to N will yield 
a local super- Hamiltonian constraint Ck + Cv = SttGJ 1 . 

Variation with respect to the shift, N' 1 , yields the su- 
permomentum constraint, 



VjTt^ = 8ttGJ 1 



(2.10) 

where the supermomentum, it ij , and matter current, J % , 
are 

5C K 



7T 



J' 



-N 



-K*' + 0—Z) Kg". 



SC M 



(2.11) 



Varying with respect to gij, on the other hand, leads to 
the dynamical equations, 

(V9^) ' = -2 (K 2 f +2(1-0 KK» 
+ i-v fc (N k Ti l 3 - ATV ife - N^ lk ) 



-C K g tJ + 



(2.12) 



where (AT 2 ) = K a Kf, and 



F" : -4 ^ y® Cv) =j29sC°(F s f. (2.13) 



y/9 Sg t 



s=0 



The constants are given by go = 2A£ 2 , g\ — —1, and 
n s = (2,0,-2,-2,-4,-4,-4,-4,-4). The stress 3- 
tcnsor is defined as 



2 S(^gC M ) 



(2.14) 



V9 S 9ij 

and the geometric 3-tensors (Fg)^ are defined as follows: 



(Fa) i:i 

(F2) ij 
(F^j 

(F 5 ) t1 



- 2 9ij, 



R 



V 2 > 



2 (Rij - ViVj) R - - gij (R - 4V 2 ) R, 
V 2 Rij - (ViVj -3R ij )R- i(R 2 ).. 

+^9ij(SR k iR kl + y 2 R-2R 2 ), 

3 (Rij - V,-Vj) R 2 - ^ gij (R - 6V 2 ) R 2 , 
{Rij + ViV,-) (RkiR kl ) + 2R (R 2 )^ 

+V 2 (RR^) - V k [V, (RR jk ) + V, (RRik) 
-~ 9ij [(R-2V 2 ) (R k iR kl ) 
-2V k Vi (RR kl )] , 



{F,\ 



3 " 

ij + 2 . 

V* 1 (V ? - (i? 2 



3(i? 3 ) 



V 2 (i? 2 ) 



j k 



1 

- 2 9ii 



I 1 L m k 

2 i 



i, \ / 



R k RL,RT — 3Vi.V; (i? 2 ) 



(J?V)„ = 2ViVj (V 2 i?) - 2 (V 2 i?) R 



+ (V<i2) (Vj-iJ) - ~<fe [(Vi?) 2 + 4V 4 i? 
(F 8 ) v . - V 4 ^- - V fc (V.V 2 ^ + V, V 2 ^) 

- (ViRf) (VjRl) - 2 (\7 k R\) (VkRjt) 



;9ij 



(V k R h 



2 (V fc V,V 2 i? fen 



(2.15) 
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The matter quantities (J*, J\ r lJ ) satisfy the conser- 
vation laws @, [23, 



d 3 x^g 



9kir 



1 

71 



2N k . , • 



0, (2.16) 



(V fe 7V, - V,7V fe ) = 0. (2.17) 



It should be noted that the energy-momentum tensor 
in GR is defined as, 



1 



6g® 



(2.18) 



where fi,v = 0,1,2,3, and g$ = -N 2 + N*Ni, g$ = 

Ni, and g^j — gij. Introducing the normal vector n v to 
the hypersurface t = const as, 



n. 



one can decompose as [29j, 

Si 

s 



(2.19) 



T tlI/ n tJ, n u , 



(2.20) 



where ft^D is the projection operator, defined as /i^ 4 J = 
.9^ + n^n u . In the GR limit, one may identify J*, Jj, Ty- 
with —2pn, Si, Sij, respectively. 



III. SPHERICALLY SYMMETRIC STATIC 
SPACETIMES 

The general spherically symmetric spacetime that pre- 
serves the form of Eq. (|2.1[) with the projectability con- 
dition is described by the metric, 

ds 2 = -N 2 (t)dt 2 + e 2l/ ^ r \dr + N r (t,r)dt) 2 

+R 2 {t,r)dVt 2 , (3.1) 

in the spherical coordinates x l = (r,6,(j>), where dVL 2 
'ill' + sin 2 9dcj> 2 and N* = {N r ,0, 0}. Clearly, it is invari- 
ant under the transformations, 



t = f(t% r = g(t',r'), 



(3.2) 



where / and g are arbitrary functions of their indicated 
arguments. With this gauge freedom, we see that, with- 
out loss of generality, we can set 



N(t) = l, R(t,r)=r, 



(3.3) 



a gauge we refer to as the canonical ADM gauge. From 
now on we shall work with this gauge. 

To consider spherically symmetric static spacetimes in 
the HL theory with projectability, we assume that there 
exists a timelike Killing vector, along t, namely £ = 
dt- It can then be shown that the Killing equations, 
Zn;v + = 0, lead to 

u(t,r)=v(r), N r (t,r) = N r (r), (3.4) 

for which the metric can be finally written as 

ds 2 = -dt 2 + (efWdt + e v ^dr) 2 + rW, (3.5) 

where 

H = v + to.N r , N r = e"- U . (3.6) 
For the metric (13.51). we find 







R i:j 


2z/ 

crrr 
= —°i°3 


+ e 


C K 


= g 2(A»-f) 


Cm' 2 


C v 


- E4 S) 

s=0 





4(1-0 1 2(1 - 2fl 

t 1 o 



(3.7) 



where a prime denotes the ordinary derivative with re- 
spect to its indicated argument, £2y = 5f5® + sin 2 OSfsj, 

and Cy s are given by Eq. (|A.l[) . Then, we find that the 



Hamiltonian constraint 



reads. 



J (C K + C v ~ 87rGJ*) e v r 2 dr = 0, 
while the momentum constraint (|2.10[) yields, 



(3.8) 



^' + (/x'-i/y 



2£ 



where 



(3.9) 
(3.10) 



J 1 = e-^ +u \v,0,0). 
The dynamical equations (I2.12[) . on the other hand, yield, 

e(2//'-2^y + Ai ' 2 ) 

4 r / / . x ,1 2(1 - 4£) 

= -2e 2(l/ -^ (V 2l/ F rr + 87rGp r V(3.11) 
2(1-0(^"-mV) + (4-3£K 2 



2(1 - 2C) 



(2,/-*/) 



-2e 2 ^ ^F gg + SnGpe^j , (3.12) 
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where 



(3.13) 



and Fij is given by Eqs. (|2.13[) and (|A.2[) . In this paper 
we define a fluid with p r = pg as a perfect fluid, which 
in general conducts heat flow along the radial direction 

Since the spacetime is static, one can see that now 
the energy conservation law (|2.16l) is satisfied identically, 
while the momentum conservation (|2.17|) yields, 

2 

Vfl' - (V - p' r ) (y - p r + p s ) = 0. (3.14) 



where pn, q, p r and pg denote, respectively, the en- 
ergy density, heat flow along radial direction, radial, and 
tangential pressures, measured by the observer with the 
four- velocity n M . Then, combining it with Eq. (|2.20l) one 
can see that such a decomposition is consistent with the 
quantities J', J 1 and T\j defined above with v = ge M . It 
should be noted that the definitions of the energy den- 
sity ph, the radial pressure p r and the heat flow q are 
different from the ones (p a , pn, q ) given by Eq. (|B.7[) . 
which are defined by assuming that the fluid is comoving 
with respect to the orthogonal frame (|B.5|) . 



It should be noted that this equation is not independent 
of Eqs. (13J1) . (|3.11l) and (|3.12l) . As a result, one cannot 
use it as an additional condition to determinate the six 
unknown functions, p, J*, v, p r , pg. But, since it 
involves only first-order derivatives, it is often very use- 
ful to use it to replace one of the three equations, (|3.9[) . 
(13. lip and (|3.12p . In summary, in the present case we 
have one integral Hamiltonian constraint (|3.8[) . one mo- 
mentum constraint (|3.9p . and two dynamical equations, 
(13. lip and (I3.12p . for the six unknowns. Thus, in order 
to determinate them uniquely, additional conditions are 
required. 

In GR, for a perfect fluid [cf . Appendix B] , one condi- 
tion usually comes from the equation of state, p — p{p), 
where p denotes the energy density of the fluid and is 
related to J 1 (but not exactly equal to it Q)- However, 
in the present case, the Hamiltonian constraint is not lo- 
cal and cannot be used to close the system. Therefore, 
to have the system closed one may take the Hamiltonian 
constraint (|3.8p as a constraint on J*(r) and then use 
the three remaining differential equations (|3.9p . (13. 
and (|3.12p . together with two additional conditions, to 
determinate uniquely the five unknowns, p, z/, v, p r and 
pg. In this paper, we shall adopt this strategy and obtain 
one of the two conditions by specifying the spatial cur- 
vature. In particular, we will consider two cases: (i) the 
spacetime is spatially Ricci-flat, and (ii) the spatial cur- 
vature is a nonzero constant. Certainly, one can equally 
choose other physical conditions to close the system. 

To relate the decomposition of the quantities J', J 1 
and Tij defined above to the ones introduced in Eq. (12.201) , 
one may introduce another spacelike unit vector, Xni 
which is orthogonal to n^, 8^ and </>^, where 



(3.15) 



(3.16) 



Such a Xfj, is uniquely determined as 

In terms of the four unit vectors, n M , x M , M , and 4>^, 
one can decompose the energy-momentum tensor for an 
anisotropic fluid with heat flow as 



+PrXt*Xv + Pe {6y,0v + <t>u.<t>v) , (3.17) 



IV. SPATIALLY RICCI-FLAT SOLUTIONS 



Requiring that the spacetime be spatially Ricci flat, 



Rij — 0, we find that v — 0, and 



<„' 2 _ 4 ( X ~ ,/ _ 2 (! ~ gO 



(4.1) 



Then, Eqs. (fA~T|) and (fA~2|) yield 



C v = 2A, Fa = ~kg l3 . (4.2) 



Inserting the above into Eqs. ipEgj) . (j3~§j) . (j3~TT|) . (|3~T2)) 
and (|3.14p . we obtain, respectively, 

J (C K + lQnGp)r 2 dr = 0, (4.3) 

t(p" + p' 2 ) + ^p! - § = -8^Ge- 2 ^, (4.4) 

e(2 M " + M '- 



." , , 4 ../ , 2(1 -4Q 

r 2 



= -167rGe- 2 "(p P +PA), (4.5) 
2(l-£)M"+(4-3^' 2 + 4 ii-^' 

= -W7rGe- 2 ^(pg+p A ), (4.6) 
2 

Vfl' - (V - p' r ) - -(v-p r +pe) = 0, (4.7) 



where 



J 1 



A 



p = PA --, PA = - PA = — 



(4.8) 



To study the above equations further, we consider the 
cases £ = and £ ^ separately. 



A. 5 = 



When £ = 0, from Ea.(|4~4|) we find that v = 0. Then, 
from Eq. (|4.ip we obtain 



£ K = -2(2r// + 1) 



a 2n 



(4.9) 
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while Eqs. (|4.5[) - (|4.T[) reduce, respectively, to 

2r M ' + 1 = -8irGr 2 e- 2f * (p P + p A ) , 
2 , 



(4.10) 



//' + 2// 2 + -fjf = -MGe-^ipe +p A ), (4.11) 



Pr + - (Pr-Pe) = 0. 



1. de Sitter Schwarzs child Solution 



If we further require the fluid be perfect, p r = Pe = Pi 
Eq. (|4.12[) tells us that the pressure must be constant. 
Without loss of generality, we can absorb this constant 
into pa, and then Eqs. (I4.5[) and (|4.6p have the solution, 



1, ( M A r, 

/i= 2 ln V + 3 r 



(4.13) 



Inserting it into Eq. (|4.1j) we find that = — 2A. Then, 
Eq. (|4.3[) requires J* = 0. This is exactly the de Sitter 
Schwarzschild solution written in the ADM form. It was 
first noticed in the framework of the HL theory in [l5j . 
and rederived later by several others. When A < 0, it 
represents the anti-de Sitter Schwarzschild solution. 



2. Anisotropic Fluid 

When pe = 7p r , where 7 is a constant but 7 ^ 1, from 
Eq. (|4T2]) we find that 



(4.14) 



where cq is an integration constant. Substituting it into 
Eqs. (|47TU| and (|4TTj) . for A = we obtain 

2?y + l = ScoTrGe-^r 2 '' , (4.15) 
r 2 fi" + 2r 2 fi' 2 + 2rp' = -Sco^Ge"^^ 7 , (4.16) 



from which we find that 



/' 

Mo 



1 



M /r\ 



In 



r\ 2 t 

l + 2 7 
' 8c 7rG£ 2 - 



(4.17) 



Clearly, when 7=1, choosing Co = — A/(87rG) and I = 
•y/3/A, the above solutions reduce exactly to the (anti-) 
de Sitter Schwarzschild solution (14.13[) . When 7^1, for 
the solutions to be real, we must have 



for £ 2 t > 0, and 



< -i, c > 0, 
> -§, c < 0, 



> -i, c > 0, 
< -|, c < 0, 



(4.18) 



(4.19) 



for £ 2 ' f < 0. On the other hand, inserting the above 
solution into Eq. (|4.9|) we obtain 



C 



K 



167rGco 

r 2(l-7) ' 



(4.12) Setting 



where 



p(r) = -cor 2 ^" 1 ) + p(r), 



^00^°°, r — > 00, 



(4.20) 



(4.21) 



(4.22) 



with p Cl do, Oqo, /?o and /3oo being constants, we find 
that the Hamiltonian constraint (|4.3p requires 



(4.23) 



/So > -3, Poo < -3. 



However, to have the center, r = 0, free of spacetime 
singularity, from Eqs. (|4.14l) . (|4.21[) and (|4.22l) . we find 
that we must assume that 



7>1, /?o>0. 



(4.24) 



Eq. (|4.18l) then shows that this is possible only when 
co < if l 2 " 1 > 0. Hence, the pressures become negative 
in this case. However, unlike GR [cf. Eq. (|B.25[) ]. all 
the three energy conditions: weak, strong and dominant 
30], can be satisfied by properly choosing p(r), which 
now is only subject to the global Hamiltonian constraint 
(|4.3p . On the other hand, when i 2 " 1 < 0, the condition 
(I4.24[) can be satisfied for Co > 0. Then, both the energy 
density and pressures can be positive. Once again, due 
to the integral form of the Hamiltonian constraint, one 
can always choose p properly, so that all the three energy 
conditions can be satisfied. It should be noted that the 
above conclusions do not contradict with the results ob- 
tained in [25[ , in which perfect fluid with some conditions 
between p and p was considered. 

When r — > 00, the pressures become infinitely large, 
and a spacetime singularity is indicated to exist there. 
One may cut the spacetime at a finite radius, and then 
join the solution to the de Sitter Schwarzschild solution, 
given by Eq. (|4.13[) . We shall consider this issue in Sec. 
VI. 



b. e^o 

When £ 7^ 0, to have the system (|4.4j) - (|4.6p closed, one 
additional condition is required. In this paper, we shall 
take it to be v = 0. Certainly, other physical conditions 
might be equally possible. Setting v — 0, we find that 
Eq. (14. 4p has the general solution, 



fi = In [ar + — 



(4.25) 
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where a and b are integration constants. Then, for A = 
Eqs. (|4~5j) and (|4l)|) yield, 



16ttG 
3 

16ttG 



a 2 (3£-2) 



r 

^•6 



(4.26) 



while Eq. (|4.7|) is satisfied identically, as one would ex- 
pect. On the other hand, for the solution (|4.25p . we find 
that 



Thus, setting 



3a 2 

P = Tl^+PV')' 



inG 



(4.27) 



(4.28) 



we find that the Hamiltonian constraint (|4.3p reads, 

p(r)r 2 dr = 0. (4.29) 



From Eqs. (|4.26l) and (|4.27l) we can see that the space- 
time is usually singular at the center, r = 0, unless 6 = 0. 
In the latter case, we have 



Pr = Pe 



3(3$ -2) a 2 f > 0, £ > 2/3 
16ttG 1 < 0, £ < 2/3. 



(4.30) 



That is, now the fluid is a perfect fluid with constant 
pressures, which are non- negative for £ > 2/3. Setting 
p(r) = 0, for which the Hamiltonian constraint is satisfied 
identically, we find that the energy density also becomes 
a positive constant. 

It is interesting to note that p r = pe = when £ = 2/3. 
In other words, the fluid becomes a dust. In GR, a dust 
cannot have a static configuration [26?] . But, now the field 
equations involve second-order derivatives of p (in GR 
only the first-order terms are involved.), which produce 
a repulsive force to hold the collapse [3l|. If we further 
choose p = — So 2 / (%ttG) so that p = 0, the corresponding 
spacetime becomes vacuum. 



V. SOLUTIONS WITH NON-ZERO CONSTANT 
CURVATURE 



which has a particular solution v = 0, k = 0. These are 
the solutions studied in the last section. Therefore, in the 
following we shall consider only the case where v' ^ 0. 
Then, we find that 



1, / k , 

-In 1 r 2 

2 V 6 



(5.3) 



for which we have 



Rij — 
C v = 



C 



2A-fc+ 3g2 + g3 fc 2 



3C 2 

9.94 + 3^5 + 56 , 3 



K 

where 
Jo = -A + '- 



-r 2 e 2 » 



9C 4 

Hp' 2 



4(1-0., 



-/' 



2(1-20 



k 3g 2 + 93 , 2 , 9ff 4 + 3g 5 + g 6 u3 
6 18C 2 18C 4 



(5.4) 



fc 3 . (5.5) 



To study this case further, we again consider solutions 
with £ = and £ 7^ separately. 



A. £ = 

When £ = 0, we find that the corresponding dynamical 
equations, momentum constraint and the conservation 
law can be written, respectively, as 



'l_fc\ , fl_ k 
r 6 / \r 2 6 



= -e" 2 " ( J + ZvGpg) , (5.7) 

(5.8) 



-2e-^ (Jo + 87rGp r ) , (5.6) 



24ttG 



e 2 ", 



Pr + \ (Pr-Pe) = ~^- r (rp' + 2) • (5.9) 



From Eq. (|3~7)) . we find that 

2e- 2 " 



R = 



2ru' - (1 - e 2v ) 



(5.1) 



When R is a constant, say, fc, the above equation can be 
cast in the form, 



2rz/ 



1 - y ) e - 1 = °' 



(5.2) 



1. Vacuum Solutions 

When the spacetime is vacuum, v = p r = pt 
the above equations show that we must have 



/1 = —00, Jo = 0, 



for which we have 



N r 



0, £k = 0, 



J* = 0, 
(5.10) 
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A 



352 + 33 , 2 , 9 34 + 3g 5 + g 6 



18C 2 

2fc , 4(352 + 53) ,2 

T + 9C 5 



18C 4 



2(9g 4 + 3g 5 +g 6 ) j3 
9C 1 ' 



(5.11) 



The Hamiltonian constraint (|3.8[) will be satisfied identi- 
cally if the coupling constants are chosen so that Cy = 0. 
These solutions are (Einstein) static universe solutions 
for either sign of the spatial curvature k [24| • 



2. Perfect Fluid 



For a perfect fluid, Eqs. (|5.6[) and (|5.T[) yield 



1 *r»W + 2 M *)-£M'-3-£= - (5-12) 



Setting 



^r]+ilnw(z), z = Jl-^r 2 , (5.13) 



we find that Eq. (|5.12[) can be cast in the form, 



(l-z 2 )w" -2zw' + 
with 



»(o+l) -- 



w = 0, (5.14) 



1 „ , 3 
- + 2i, b= -. 

2 ' 2 



(5.15) 

The general solution of Eq. (|5.14[) is given by 

w = Cl P b a {z) + c 2 Q b a (z), (5.16) 

where c\ and ci are the integration constants and should 
be chosen so that the solution is real. P b {z) and Q a {z) 
are, respectively, the associated Legendre functions of the 
first and second kinds [33]. Inserting Eq. (|5.16l) into Eq. 
([5TT3|) . we find that 

H = \ ln(r) + ~ In [ Cl P b (z) + CaQfc)] + Mo , (5.17) 

where fi = [ln(fe/6)]/8. Then, Eqs. ([377]) . (l5To) and ([575]) 
yield, 



A" 



3( Cl p a b (z) + C2 g^(z) 



3A/l-|r2 



Cl P a b (z) + c 2 Q a b (z) 



8ttG 



,2^0 



16ttG 



^+6M C l P aW + C 2Q b W, 



^/i-^( Cl p; b (z)+ C2 Q; b (z) 



(5.18) 



Note that in writing the above expressions we did not 
impose any conditions obtained from the vacuum case, 
so that our solutions are as much applicable as possible. 
In particular, by properly choosing the parameters, we 
can have To < so the pressure in the center of the star 
is positive and the fluid satisfies all the energy conditions 
[See the discussions below.]. On the other hand, setting 

p{r)=p k +p{r) 



16ttG 

we find that the Hamiltonian constraint 
p(r)r 2 dr 



1 



6 r2 



0, 



where 



Pk 



16ttG 



, _ 3ff2 + 93 , 9.94 + 3.95 + . 



3C 2 



9C 4 



(5.19) 
reads 

(5.20) 



(5.21) 



To study the singular behavior of the solution near the 
center, we first note that [Hf 



1 



r(i-6) \z-i 



1 



6/2 



x F ( —a, a + 1:1 — b: 



b/2 



Z- 1 



^(6)^^-0,0+151-6;^ 



r(l + a + 6)r(-6) fz - 1 



2T(l + a-6) 
xF ( -a,a+ 1; 1 + b; 



z + 1 

1 - z 



, (5.22) 



for |1 — z\ < 2, where F(a,b; c; z) denotes the hypergeo- 
metric function. When \z\ <C 1, it is given by 

t-,/ 7 ^ -. a b ab(a + l)(b+l) , „ , , x 

F(a, b; c; z) ~ 1 + — z + v . -z 2 + O (z 3 ) . 

c c(c+l) v 7 

(5.23) 

Thus, as r — > + , we find 



p*( . gj_jl 17fc r 2 425fc2 

- r 3/2^ 48 4 6 08 



O 



«2 
r 3/2 ' 



17k 
48~' 



+ a 3 r 3 + o(r 5/2 ) 1,(5.24) 
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where 



a i 



0,3 



2 3 / 2 V2T (| 

a 2 = — i 



fc\3/4' 



r(-i) (-1) 
r(-f)r(2 + 2i) 
r(§)r(-i + 2i) V 24 



k\3/4 : 
6 I 



3/2 



(5.25) 



Then, from Eqs. (j5"T8|) and (l5T9|) we find that when 

a\C\ + a 2 c 2 = 0, (5.26) 

all the quantities, Ck, p, v and p are free of singularity 
at the center r = 0. In fact, for such a choice we have 



C 



K 



-6c 2 a 2 a 3 e 2f *° , t> ~ 0, 



87rG" 



(5.27) 



r -> 0+. 



b. e^o 

When £ ^ 0, the dynamical equations and the momen- 
tum constraint read, respectively, 



(V + m' 2 



4 (2 + £)fc 



r 3 
2(1-40 , (1 + 20* 



1 - -r" 



= -2e~ 2 '' (Jb + 8irGp r ) , (5.28) 
2(1 -£)//' + (4 -3£)m' 2 

4(1 - 2Q (3 - 5£)fc 1 , (1 - 2£)fc 

r h 

r 3 J P 3 

= -2e- 2 " (J- + 8irG Pe ) , (5.29) 

e(i-W (m"+m' 2 )+£(^> 



2£ , fc 



-87rGe- 2 ^w. (5.30) 



by 



1. Vacuum Solutions 



In the vacuum case, one of the solutions is still given 



/i = — oo, J-q = 0, 



(5.31) 



for which we have N r = = Ck- When p, is finite, from 
Eqs. (|5.28[) - (|5.30|) we find that the vacuum equations 
can be cast in the forms, 



4 (2 + 0* \ , 



e(l-=r»j(2 M " + M*) + (=---3 — '"J/' 



2(1-40 , (1 + 20* 



= -2er 2 »F 0) (5.32) 



1 - -r' 



and 



2(1 - £),/' + (4 - 30m' 2 
4(1 - 2Q (3 - 5£)fc 1 , 

q r M - 

r 3 
1,. _ ,x 3 2(2-0 



(1-x 2 )m,xx + ~(4-5:e 2 )m,x j - 



(1 - 2Q_fc 
3 

(5.33) 
= 0, (5.34) 



where x = \/k/6 r. Eq. (|5.34l) has the general solution 
3(3£-2) 



Mo + ln(r) 



Ml 



6 ' 



2£/b 
2-3^ 
""4T ' 



-arcsm 



— r 
6 , 



(5.35) 



where /^o and Mi are integration constants. Inserting the 
above into Eq. (|5.32p . we find that it is satisfied only 
when 



Mi = 0, £ 



and for which Eq. (|5.33p gives JF = 0. That is 

a * 3g 2 + gg , 2 9ff 4 + 3g 5 + g 6 3 _ 
A 6 ^SC 2- 18C 3 " 



(5.36) 



(5.37) 



When Ml = and £ = 2/3 it can be also shown that 
Ck = 0. The Hamiltonian constraint is then satisfied 
identically, when Cy = 0, i.e., 

2A-k+ 3g2 ± 93 k 2 + 9g4 ± 3g5 ± 9e e = 0, (5.38) 



3( 2 



9C 4 



as one can see from Eq. (|5.4p . Therefore, provided that 
the coupling constants g n , (n = 2, 3, 6) are chosen so 
that Eqs. (|5T37| and (fOB"]) hold, the solution 



In 1 



1 

~2"~V~ 6 
ln(r)+Mo, (£ = 2/3), 



(5.39) 



represents the unique vacuum solution of the HL theory 
with maximal symmetry for any given curvature k and 
nonzero £. It is interesting to note that £ = 2/3 is the 
case where an anisotropic Weyl symmetry exists in the 
UV limit [l|. Note also that the spacetime described by 
the solution (15.39[) is not flat even in the sense of the 4- 
dimensional geometry. For example, the corresponding 
4-dimcnsional Ricci scalar is given by, 

R (i) = 12e 2 M0 L _ + fcj ^ = 2/3), (5.40) 

which shows that the spacetime is not flat even when 
k = 0. 
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2. Perfect Fluid 



On the other hand, for a perfect fluid Eqs. (|5.28l) and 
fl5H| yield, 



k , 

1 - -r 2 
6 



(1 - 20a*" + 2(1 - 0/' 



4£ | (l~6Qfc ] 
r 6 



3=0. (5.41) 



Setting 



1 + 2g 
4(1-0 



ln(r) 



1-2^ 
2(1-0 



hiw(r), 



(5.42) 



we find that Eq. ()5.41|) can be cast in the form of Eq. 
(|5.14p . but now with 



1-2$ + + 
2(2£-l) 



»4 



(5.43) 



Therefore, in the present case the general solution of 
(|5.41[) is given by 



l + 2£ 

4(1-0 



ln(r) 



1 - 2g 
2(1-0 



where, as previously, z 



In [c 1 P b a (z) + c 2 Q a b (z)\, 

(5.44) 

1 — |r 2 . Once fi is given, 

from Eqs. (|5.28[) and (|5.30p we can find p and v, which 
are too complicated to be written explicitly here. 

To study the asymptotical behavior of the above so- 
lutions near the center, we first notice that Pa{z) and 
Qa( z ) take the same forms as those given by Eq. (|5.24l) . 
as they do not depend explicitly on the parameter a as 
r -> 0. We find that 



Jo , 301-40 (aid 



5 



a 2 C 2 



8ttG 64ttG(1 
901 - 20 (aid a-u-2 



32ttG(1 - 2 



1-5 



c 



K 



305 - 80 ( aid + a 2 c 2 
4(1 -0 2 



(5.45) 



as r — >• 0. Thus, when —1/2 < £ < 1, all these quantities 
are finite at the center for any given d and c 2 , provided 
that aid + a 2 c 2 ^ 0. When £ > 1 or £ < 1/2, they 
diverge there unless d and c 2 are chosen such that aid + 
a 2 c 2 = 0. Therefore, in the present case d and c 2 must 
be chosen so that 



aid + a 2 c 2 



#0, -1/2 < f < 1, 
= 0, otherwise. 



(5.46) 



VI. JUNCTION CONDITIONS 

Let us consider a surface E, defined by r — r , in the 
spacetime described by the metric (I3.5[) . which divides 



the whole spacetime into two regions, the internal region 
r < r Ql and the external region r > r , denoted, respec- 
tively, by V~ and V + . Note that once the metric is cast 
in the form (|3. 51) , the coordinates i and r are all uniquely 
defined. As a result, the coordinates used in V + and V~ 
must be the same, i.e., 



(6.1) 



Since the highest order of derivatives in the HL theory 
is six, one may require that the metric coefficients be at 
least C 6 ; that is, their derivatives up to six-order exist 
and are continuous across E. However, this requirement 
often is too strict, and, in particular, will exclude the exis- 
tence of infinitely thin shells [33J . To relax this condition, 



from Eqs. (I2.6[) and (|2.15l) we can see that the quadratic 
terms of the highest derivatives are third-order, so we 
may require that the metric coefficients be at least C 3 . 
In this way we can avoid terms that are powers of Dirac 
(5-functions, which mathematically are not well defined. 

For the spherically static spacetime, this condition is 
still very strict, since the quadratic terms of the high- 
est derivatives now are only terms involving v" 2 , i/'i/", 
and // 2 , as we can see from Eqs. (I3T71) - (|3~T2"1) . (|A.1[) 
and (IA.2|) . Therefore, without loss of generality, we shall 
assume that v(r) and /i(r) are at least C 1 and C°, re- 
spectively, across the surface r = tq, and at least C A and 
C 1 elsewhere. Denoting quantities defined in V + (V~) 
by F + (F~), we find that /i and v can be written as 



Fir) = F + {r)H (x) + F~(r) [1 - H (x)] , 



(6.2) 



where F = (//, v), x = r — ro (It must noted that there is 
no confusion between x used in this section and the one 
used in Sees. IV and V.), 

limit r ^ r +Ai + (r) = limit r ^ ? ,-^"(r), 
limit r ^ r +v + {r) = limit r ^ r -v~ (r), 
limit r ^ r +i/+(r) = )xaAt r _^ r -v~{r), (6.3) 

and H(x) denotes the Heavside function, defined as 
which has the properties , 



H n (x)=H(x), [1-H(x)] n = [l-H(x)} ) 
H(x) [1 - H (x)] = 0, H'(x)=6(x), (6.5) 

in the sense of distributions, where 6(x) denotes the Dirac 
delta function. Although the high-order derivatives of \x 
and v are not continuous across the hypersurface r = ro , 
we assume that they all exist and are finite in the limits 
r — > r^ 1 . Then, we find that 

At' = A* °, m" = A* D + [mT 5 ( x ) i 
v' = v D , v" = v' D , 
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V 



V 

(4) 



v"' D + {v"T S (x), 

v^ D H + [v"'Y 5 [x) + [v"Y 5' {x) . 



y (5) _ „(5)D 



8(x) 



S'(x) 



+ [v"\-5" (x) 



(6.6) 



where 

» 



= limit r _ s . J ,+ v + ^ (r) - limit r 
p(n)D = F + ^ H + F~ ( n \l — H) . 



{n \r), 
(6.7) 



Inserting the above expressions into Eqs. (PT7)) and (1A~T1) . 
we find that 



-A" 

(i) 
v 

(3) 



L, K , 



(1) D 
V ' 
(3) D 



(0) 
V 

C 
C 



C 



(0) D 



(2) 
V 



(3) Ira 



V 



— 



8g 7 e 



-61/ 



4„3 



C 4 r 



X 1/ 



(2)n 

2rV - (1 - e 2l/ ) 
"-<*(*), 



while from Eq. (|A.2[) we find that (-Fn)y's are given 
by Eq. (IA.3I) . The superindex "Im" represents the im- 
pulsive part of the quantity considered, which is usually 
proportional to 5(x) and its derivatives [cf. Eqs. (|6.9[) and 
(|A.5[) .]. Separating the nondistributional from the distri- 
butional parts of the matter content as 



J 1 


= J t,D 


+ J tJm , 




V 


= v D + 


v Im , 




Pr 


= p? + 


pl m , 




Pe 


- pU 


P'e m , 


(6.9) 



we find that the Hamiltonian constraint (13.81) now reads 



(8irGJ tJm - Cf /m ) e v r 2 dr. (6.10) 



While the momentum constraint (|3.9p and the dynamical 
equations (|3.11[) and (13. 12)) remain the same in regions 
V + and V~ , but across the thin shell at r = ro, they 
read 

t\p!]~5(x) = -8TrGe 2 ^'^v Im , (6.11) 
£[f/]~*(aO + e- 2 *F?™ = SnGe^-rip 1 ™, (6.12) 

= -8TxGe 2{v '^ V \ m , (6-13) 

where i 7 "/™ and Fq™ are given by Eq. (|A.4|) . Assuming 
that the matter content has distributional contributions 
no more singular than a (5-function, we see from above 
that in order to cancel the 5-function derivative terms in 



Ff.™ and Fg™, it is sufficient that there is some tuning 
of the couplings as <?g = 857/3. 

It is interesting to note that in the GR limits: £ = 
and C -> 00, we have F/ r TO = F*™ = 0, and Eqs. - 
(I6.13P reduce to 



Tm 



,1m 



e 2(M-") 

8ttG 

^ = 0, (£ = 0, C^OO). 



(6.14) 
(6.15) 



Thus, in this limit the radial pressure of the infinitely 
thin shell always vanishes. This is consistent with the 
conclusion obtained early by Santos [26j. 

However, this is no longer true when £ 7^ even in the 
low energy limit where i 7 "/™ = Fg™ — 0, as can be seen 
from Eqs. (|6.1ip - (|6.13p . In particular, when v = 0, we 
find that 



(1- 



-e 2f +S(x 
-e 2fl 5(x 



SttGv 1 " 1 , 
SirGp 1 ™, 
2 ^5{x) = -SnGp^. 



(6.16) 
(6.17) 
(6.18) 



This completes the general description of the junctions 
of a spherically symmetric star in the HL theory of grav- 
ity. In the rest of this section, we shall apply the above 
general formulas to the solutions found in the last sec- 
tions. We first notice that solutions with nonzero con- 
stant curvature k cannot be matched with the ones with 
zero constant curvature. This is because in the former 
the function v cannot be zero for any given ro. As a 
result, v cannot be continuous across r — ro- There- 
fore, only the solutions with the same curvature k can 
be matched to each other. However, since £ is a running 
coupling constant, in principle £ can have different values 
at different energies. In particular, the spacetime deep 
inside a very massive star is expected to have a very high 
temperature, and one would expect that £ will have dif- 
ferent values in the regions inside and outside of the star. 
Thus, in the following we shall consider the possibility of 
matching a fluid to a vacuum solution with different £. 
In addition, we shall consider only the match without an 
infinitely thin shell at r = ro, that is, wc shall set 



P 



Im 



,Im 



p 



Im 



0, 



(6.19) 



which implies that p, and v must be at least C 1 and C , 
respectively. 



A. Spatially Ricci Flat Solutions 

When the spacetime is spatially Ricci flat, in Sec. 
IV we showed that the de Sitter Schwarzschild solution 
(|4.13p is the unique vacuum solution. Therefore, in this 
case the spacetime outside the star is uniquely described 
by this solution, 



1 



= - In h — 



M 4 



A 



v+ = 0. 



(6.20) 
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Inside the star, two solutions were found, one is for £ = 
given by Eq. (|4.17p and the other is for £ 7^ given by 
Eq. (|4.25j) with b = 0. In the following let us consider 
them separately. 



1. £ = 

When £ = 0, the spacetime inside the star is described 
by Eq. (|4.17|) . which now can be written as 





'M_ 






h 




2 


r 





Mo, 



v- = 0. 



(6.21) 



From the above expressions we can see that v is analyt- 
ical across r = ro, while the condition that /x being C 1 
requires 



M+M = (J>-(ro), 



(6.22) 
(6.23) 



Inserting Eqs. (|6.20j) and (|6.21j) into the above condi- 
tions, we find that 



3M 4 



Ar 3 



2"; 



= 3e 2w \ 

_ M_-2 7 r (?) 27 

3M + + Ar 3 M _ +rQ (n>) 2 ^ : 

from which we obtain 



3M+ - 2Ar 



M 4 



3M. 



+ 2(l- 7 )ro (y) 



2- 



A 



SirGcor^-^ . 



3.24) 
3.25) 

3.26) 
3.27) 



Note that the condition of Eq. (|6.27p guarantees that the 
radial pressure is continuous across the surface r = r , 
i.e., p r (r ) — Pa, as can be seen from Eq. (|4.14l) . 



2. £7^0 

When £ 7^ 0, in Sec. IV, we found the perfect fluid 
solution given by Eq. (|4.25[) with b = 0, that is, 



/!_ = ln(ar), V- = 0. 
The corresponding pressure is given by 
3(3£- 2)a 2 



P = 



16ttG 



(6.28) 



(6.29) 



It is interesting to note that this solution is exactly the 
de Sitter solution in GR. However, in the HL theory it 
corresponds to a perfect fluid with its pressure given by 
Eq. (|6.29|) . As shown explicitly in Sec. IV, choosing 
p(r) = the energy density becomes p = 3a 2 /(8nG) [cf. 
Eq. (|4.28[) ] , which satisfies all the three energy conditions 



for 4/9 < £ < 4/3. For such an internal solution, the 
conditions (|6.22l) and (I6.23[) read 



Af+ = 0, 



(6.30) 



B. Stars with Non-Zero Constant Curvature 

When the spatial three-curvature R is a nonzero con- 
stant, we found two vacuum solutions, one is given by Eq. 
(I5.3ip with /1 = —00 (N r = 0), and the other is given by 
Eq. (|5.39[) with £ = 2/3. The one with (1 = —00 cannot 
be matched to any solution with finite \i across r = r . 
As a result, the only possible solution that describes the 
spacetime outside of the star in the present case is the 
one given by Eq. (I5.39[) . 



ju+ = ln(r) + MO) ^+ 



(6.31) 



On the other hand, two perfect fluid solutions were found, 
one is for £ = given by (|5.17|) . and the other is for £ 7^ 0, 
given by Eq. (I5.44j) . Redefining the integration constants 
Ci and C2 appearing in Eq. (|5.17[) . we find that in both 
cases the solutions can be written as 



l + 2£ 
4(1-0 



ln(r) 



l-2£ 
2(1-0 



In 



6 r 



C2Qi{z) 

(6.32) 



Clearly, in the present case v is analytical across r = ro, 
and we have F/ r TO = F/™ = 0. The conditions of Eq. 
(|6.19[) reduce, then, to those given by Eqs. (|6.22l) and 
(|6~23|) . For the solutions given by Eqs. (I63T1) and ([6T32]) , 
those conditions read 



Mo = 



l-2£ , Cl P*(z ) + c 2 Q b a (z Q ) 



In 



2(1-0 

Cl P' a b (z ) + c 2 Q' a b (z ) 
Cl P b (z ) + c 2 Q b a {z Q ) 



3/2 
r o 

9zo 



(6.33) 



(6.34) 



where zq = z{tq). Clearly, by properly choosing the free 
parameters involved in this model, the above equations 
will be satisfied. 



VII. CONCLUSIONS 

In this paper, we have systematically studied spher- 
ically symmetric static spacetimes filled with a fluid in 
the HL theory of gravity with projectability, but without 
detailed balance conditions. 

After writing down the relevant field equations cou- 
pled with a fluid (including the Hamiltonian, momentum 
constraints, dynamical equations, and conservation laws) 
in Sec. Ill, we systematically studied spatially Ricci flat 
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spacetimes in Sec. IV, and spacetimes with nonzero con- 
stant curvature in Sec. V, for both cases where the space- 
times are vacuum and filled with a fluid. In particular, 
in Sec. IV we showed that the de Sitter Schwarzschild 
solution is the unique vacuum solution that is spatially 
flat. In this section, we found two classes of solutions 
coupled with a fluid. The first class, given by Eq. (|4.17l) . 
represents spacetimes filled with an anisotropic fluid in 
which the tangential pressure is proportional to its ra- 
dial pressure, given by Eq. (|4.14l) . The second class is 
given by Eq. (I4.25[) with 6 = 0, which represents a per- 
fect fluid with constant pressure. This class of solutions 
actually describes the de Sitter space, but corresponds 
to a perfect fluid with positive energy density and pres- 
sure. This is in contrast to GR, where de Sitter space 
does not satisfy the strong energy condition 30]. The 
main reason is that in the HL theory the Hamiltonian 
constraint becomes a global one, and has less constraint 
on the energy density. When £ = 2/3, the pressure van- 
ishes thus representing dust. In GR, dust cannot have 
a static configuration and it necessarily develops space- 
time singularities [2t| . In the HL theory, higher order 
derivatives are present, and it is exactly the existence of 
these terms that produce repulsive forces, which prevent 
the collapse of the dust. This provides another concrete 
example in which a would-be caustic is regularized by the 
repulsive gravitational forces, created from the gradients 
of high order derivatives of curvature [3l[ . 

In Sec. V, we found that there are two different vac- 
uum solutions for spacetimes with nonzero constant cur- 
vature. One is an (Einstein) static universe, given by Eq. 
(|5~T0)) or (|5~3T]) . and the other is given by Eq. (|Q9)) . 
which has the maximal symmetry and is not flat. The 
general solutions for a perfect fluid was found explicitly, 
and are given, respectively, by Eq. (I5.17[) for £ = 0, and 
Eq. (1Q4|) for £ ^ 0. 

To construct spacetimes that represent stars, we inves- 
tigated the junction conditions across the surfaces of stars 
in Sec. VI, and obtained the general junction conditions 
with/without infinitely thin shells. It is remarkable that, 
in contrast to GR [26[ , the radial pressure of the star does 
not necessarily vanish on the surface of the star, neither 
does the radial pressure of the thin shell. This is due 
to the high order derivatives of the spacetime curvature. 
As a result, a star can be formed much easily than that 
in GR. Applying those general formulas to the solutions 
found in Sees. IV and V, we showed that this is indeed 
the case. In particular, all the internal solutions have a 
nonzero radial pressure on the surface of the star but it 
is still possible for them to be matched smoothly to a 
vacuum spacetime without the presence of an infinitely 
thin shell on the surface. 

In Appendix B, we studied anisotropic fluid with heat 
flow in general relativity with the metric written in an 
ADM form. Among other things, we showed explicitly 
that any given static solution written in an orthogonal 
form (jB.lOp can be always brought into an ADM form 
(IB.12|) with the projectability condition by the coordi- 



nate transformations Eq. (|B.11[) . However, such coor- 
dinate transformations are not allowed by the restricted 
diffcomorphisms, (|2.4[) , of the HL theory. In particular, 
the extrinsic curvature tensor Kij and the 3-dimensional 
Ricci tensor iJy are no longer tensors under the trans- 
formations given by Eq. (|B.11|) . As a result, any actions 
constructed from 2£y and Rij, such as the one given by 
(I2.5[) . are in general not invariant and therefore the trans- 
formed metric is not a solution of the projectable theory. 
In particular, a vacuum solution is no longer vacuum in 
the new frame with the projectability condition [35] ]. 

Finally, let us comment on another delicate issue, 
namely the singular behavior of the extrinsic curvature 
K for some of the solutions discussed in this paper. We 
have not used it to identify spacetime singularities in the 
present paper, because the singularities given by K seem 
not to be as serious as those given by other quantities, 
such as the Ricci curvature R, energy density p and pres- 
sure p. For example, Cai and one of the present authors 
found in [35[ that K is sing ular at r = (3M/IAI) 1 / 3 for 
the anti-de Sitter Schwarzschild solution. This singular- 
ity is absent in general relativity and the tidal forces and 
distortions felt by observers at these singularities are all 
finite. Therefore, it is not clear whethe r sp acetimes be- 
yond these points are extendable or not [36( ■ It is exactly 
due to these considerations that we did not use K to iden- 
tify spacetime singularities, although the singularities of 
K are scalar ones and cannot be removed by the foliation- 
preserving diffeomorphisms 35]. In fact, in Appendix C 
we showed explicitly for which solutions found in this pa- 
per K is regular or singular at the center. Understanding 
the nature of singularities of K is an important issue in 
the HL theory, and we wish to return to this problem in 
the near future. 
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Appendix A: Functions and (F. 



The Lagrangians C^'s in Eq. (|3.7l) for the static spher- 
ically symmetric spacetime (|3.5p are given by 
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The functions (F s ) i - defined by Eg. (|2.15l) are given by 
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where we denoted i/™) = d n v/dr n . Inserting Eq. (|6.6 
into the above expressions, we find that 
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Appendix B: Spherically symmetric and static 
spacetimes in General Relativity 

In this Appendix we will see how static spherically 



symmetric metrics in theories with unbroken diffeomor- 
phism invariance can always be brought to a projectablc 
form. Furthermore, we specialise our discussion to Gen- 
eral Relativity and see how the equations of motion read 
in the new gauge. 

The metric for spacetimes with spherical symmetry 
takes the general form, 



ds 2 = g ab dx a dx b + R 2 dn 2 , 



(B.l) 



where a,b = 0,1, and g ab and R are all functions of x° 
and x 1 , and dQ 2 = d6 2 + sin 2 6d(f) 2 . The four- velocity of 
a fluid moving radially in such a spacetime usually has 
only two nonvanishing components, 

u M = («o,«i,0 J Q), (^ = 0,1,2,3) (B.2) 
subject to the condition, 



.A 



U\U 



1. 



(B.3) 



Clearly, the metric (|B.1[) is invariant under the coordinate 
transformations, 



x = f ( x ,x J, x =g\x ,x J, 



(B.4) 



where / and g are arbitrary functions of their indicated 
arguments. Using one degree of the freedom, one usually 
sets goi — 0. When one considers a fluid, one often uses 
the other degree of freedom to choose the coordinates to 
be comoving with the fluid, so that the four-velocity of 
the fluid is given by u M oc 5®. Then, in this gauge we 
have 

ds 2 = _ e ™(r,r) dT 2 + ^r)^ + R 2 ^ ^2 > (Rg) 

with 



% = e J M . 



(B.6) 



An anisotropic fluid with heat moving along the radial 
direction takes the form, 

T^u = PoU^y-v + Prt^Tu + Pe {Q^Q V + <j>n<f>„) 

+q Q (u^rv + u^r^) , (B.7) 

where tv, 6^ and (j)^ are unit vectors, defined by 



R sin t 



(B. 



Po, PR, Pe and q a are, respectively, the energy density, 
radial pressure, tangential pressure, and heat of the fluid 
comoving in the orthonormal frame. Note that the metric 
(|B.5[) is still invariant under the rescaling, 



r = f(r'), r = g(r'), 



(B.9) 



where / and g are arbitrary functions of their indicated 
arguments. 

When the spacetime is static, <E> and R become 
functions of r only. Then, using the remaining gauge 
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freedom (|B.9[) we can always set R{r) = r, so that the 
metric finally reads 

ds 2 = -e 2 *(''W 2 + e 2 ^dr 2 + rW. (B.10) 

Let us now make the coordinate transformations, 



T = t 



J ^-2* _ i e *dr. (B.ll) 



Then, in terms of t, the above metric takes explicitly the 
canonical ADM form with the projectability condition, 

ds 2 = -dt 2 + (e^dt + e^drY + r 2 dfl 2 , (B.12) 



with 
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Clearly, to have the coordinate transformations be real, 
we must assume 



e 2 * < 1. 



'J is often written as 137 



e 2* = 1 _ 2m ( r ) 



(B.14) 



(B.15) 



where m(r) represents the gravitational mass within the 
shell r. When m(r) > 0, the condition (|B.14[) is satisfied 
identically. 

It should be noted that the coordinate transformations 
given by Eq. (|B.llj) are not allowed by the foliation- 
preserving diffeomorphisms (|2.4j) . In particular, the ac- 
tion of Eq. (|2.5[) is not invariant, because now the ex- 
trinsic curvature and Ricci tensors K y and Rij no longer 
behave like tensors under these transformations. 

Let us note here that the definitions of the energy den- 
sity p , the radial pressure pr and the heat flow q are 
different from the ones (pu, Pr, q) given by Eq. (|3.17l) . 
which arc defined by assuming that the fluid is comoving 
with respect to the canonical ADM frame (IB.12[) . The 
relation between the two sets of quantities is the following 

PH = _ 2ft (po + e 2 "p R - 2e%), 

Pr = _ ^ (Pr + e 2 > - 2e%), 

1 = j-r^r[( 1 + e a/ *)?o-^(^ +Pfl)].(B.16) 



The nonvanishing components of the Einstein tensor 
for the metric (|B.12p are given by 
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Then, for an anisotropic fluid (IB.7|) . the Einstein field 
equations, G^ v = 8irGT^ v , yield, 

2rp! - 2r (1 - e" 2 ") v' - (l - e 2 ") e" 2 " + 1 

= 8TrGr 2 e 2 ^-^ Po , [Goo], (B.18) 

(1 - e- 2 ") 1/' = -^Gre^"-") (p c - 2e%) , 

[e 2 ^Goo + (l-e 2M )G 11 ], (B.19) 

r (//' + 2^ 2 - iV) + 2m' - (1 - e" 2M ) ^' 

= -8irGre 2 ^-^ Pei [G 22 ], (B.20) 

9o = 0, [e^Goo + (l-e 2 ^)Goi]. (B.21) 

The last equation shows clearly that in GR heat flow 
along radial direction is not allowed in static spherically 
symmetric spacetimes. 

The conservation laws V^T^ = 0, on the other hand, 
give 



2rq p! + (1 - e" 2 ") (rq' + 2q ) = 0, 



(B.22) 
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(B.23) 



For q a = 0, Eq. (|B.22I) is satisfied identically, while Eq. 
(|B~23t reduces to 



{po+PR)n'+(l-e- 2 ») [PR + -(PR-Pe)\ =0, (q o = 0). 

(B.24) 

When v = 0, from Eq. (|B.19|) we see that p# = — p Q . 
If in addition we have a perfect fluid pr — pg, from Eq. 
(|B.24I) we find that p' R — 0, that is, the pressure is con- 
stant. Then, from Eq. (|B.18|) we find that p is exactly 
given by Eq. (I4.13[) . which is identically the de Sitter 
Schwarzshcild solution written in the ADM form [Hj]. 

When v = and p g = JPr, from Eqs. (|B.18|) - (|B.24|) 
we find that 



" = 2 ln 



M 
r 



r\ 2 T 



"Pi? 



-7 Pe 



+ ^0, 

-cor 2 ^- 1 ), (B.25) 
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where fi is given by Eq. (|4.17|> . and cq is a constant. 
To have p a non-negative we must assume Cq < 0, while 
finiteness at the center requires 7 > 1. As a result, we 
find that p + pg = (1 — l)po < 0, that is, the fluid does 
not satisfies the weak energy condition [30(. In fact, it 
does not satisfies any of the three energy conditions. 

Appendix C: Singular Behavior of Extrinsic 
Curvature K 

For the static spacetimes described by metric (|3.5[) . we 
have 

K = e»~ v (V + ^ . (C.l) 

For the (anti-) de Sitter Schwarzschild solutions (|4. 1 3[) , 
K is given by [35j 

K= y^2^{ 4 - smtp )■ (c - 2) 

As noticed in [35], K is singular for the anti-de Sitter 
Schwarzschild solution (A < 0) not only at the center 
r = but also at r = (3M/IAQ 1 / 3 . The latter singularity 
is absent in GR. 

For the Ricci-flat solutions given by Eq. (|4.17[) . we 
have 




(C.3) 



which is singular at the center, unless M = and 7 > 1. 
For the solutions of Eq. (|4.25j) with b = 0, we find that 

K = 3a, (C.4) 

and therefore everywhere finite. 

For the solutions of Eq. (|5.17|) . we find that 

K = 3e M ^^(l + 0(r)), (C.5) 

and therefore finite at r = 0. 

For the solutions of Eq. (|5.39j) . we find that 



K = 3e W) y 1 - -r 2 , (C.6) 

and therefore everywhere finite. 

For the solutions of Eq. (I5.44j) . we find that as r — > 

3 , , 1-2; 3(2{-l) 

K ~ 2(1-0 1 +a2 ° 2 ^ r^zr, (C.7) 

with non-singular corrections. Thus, when a\C\ +a^C2 ^ 
and 1/2 < ^ < 1, K is finite at the center. For 
a\C\ + CL1C2 — one obains the same condition for £, but 
these solutions are forbidden from the regularity condi- 
tion (|Q6|) . 

We have not provided here the singular behavior of 
the invariant KijK lJ , but it turns out to have similar 
behavior as the one of the trace of the extrinsic curvature 
K. 
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